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Interactions between atoms and lasers provide the potential for unprecedented control of quantum states. 
Fulfilling this potential requires detailed knowledge of frequency noise in optical oscillators with state-of-the- 
art stability. We demonstrate a technique that precisely measures the noise spectrum of an ultra-stable laser 
using optical lattice-trapped 87 Sr atoms as a quantum projection noise-limited reference. We determine the 
laser noise spectrum from near DC to 100 Hz via the measured fluctuations in atomic excitation, guided by 
a simple and robust theory model. The noise spectrum yields a 26(4) mHz linewidth at a central frequency 
of 429 THz, corresponding to an optical Q of 1.6 x 10 16 . This approach improves upon optical heterodyne 
beats between two similar laser systems by providing information unique to a single laser, and complements 
the traditionally used Allan deviation which evaluates laser performance at relatively long time scales. We use 
this technique to verify the reduction of resonant noise in our ultra-stable laser via feedback from an optical 
heterodyne beat. Finally, we show that knowledge of our laser's spectrum allows us to accurately predict the 
laser-limited stability for optical atomic clocks. 



The utility of a precision frequency source is limited by its 
instabilities. For this reason, many methods to rigourously 
characterize these instabilities have been developed [1]. Ultra- 
stable lasers pose a unique challenge to characterizing fre- 
quency instabilities because, until now, measurements of their 
performance required an optical herterodyne beat between 
two or more lasers [2-7]. Single laser performance can be in- 
ferred from a three-cornered hat type measurement [8, 9], but 
valuable information about a laser's frequency noise power 
spectral density (PSD) [10] is limited in an optical beat by the 
less stable laser. 

Optical lattice trapped 87 Sr atoms are uniquely suited for 
laser noise spectral analysis due to the ultra-narrow linewidth 
and field insensitivity of the 1 S'o (\g)) to 3 Pq (|e)) clock tran- 
sition as well as the low quantum projection noise (QPN) 
one can achieve with ensembles of many atoms. To accom- 
plish this, we adopt a technique similar to radio-frequency- 
based dynamical decoupling [11-14] to manipulate the fre- 
quency noise sensitivity of the 87 Sr optical clock transition. 
We demonstrate the ability to characterize the noise spectrum 
of an ultra-stable laser with fractional frequency stability at 
1 x 1CP 16 . To guide and interpret our experimental mea- 
surements, we develop a simple and robust theoretical frame- 
work that combines concepts from the characterization of fre- 
quency instabilities [1] with a model for atomic sensitivity to 
frequency fluctuations [15-18]. We compare experimentally 
measured fluctuations in atomic population to our theory and 
accurately determine the PSD of our laser. As laser stability 
advances, we can continue to leverage the QPN-limited noise 
floor of this technique to analyze lasers with greater stability. 
At the same time, atoms could serve to stabilize lasers at the 
level of intrinsic atomic coherence [19, 20]. 

To model the frequency of our ultra-stable laser, we con- 
sider a fixed frequency with a small, time-dependent noise 
term: = uj^ a + Sui(t). The instantaneous phase of 

the laser is then given by 4>L{t) = f dt'u)z(t') — ujL t + 
f dt'5u)(t') = ujL t + 5<f)(t). The resulting Hamiltonian for a 



two level atom, with energy spacing hui a , driven by this laser 
is given by [21] 
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where A = W£ — w a , a z is a Pauli spin matrix, and fl(t) is the 
Rabi frequency. The time dependence of fl(t) relates to the 
chosen spectroscopy sequence. In the absence of other pertur- 
bations, the atom behaves as a filter, through which laser insta- 
bilities can be viewed. To calculate the output of our "atomic 
filter," we will need to know the atom's impulse response. 
This quantity describes the atomic response to frequency devi- 
ations of the form 5uj(t, to) — A(f)5(t — to), where 8(t) is the 
Dirac delta function and Acf> is infinitesimally small. Using 
the definition of 54>(t) from above, we see that this is equiv- 
alent to considering an instantaneous step in 8qb(t) at time to- 
5(f>(t) can also account for deliberate phase shifts of the laser 
which we implement during spectroscopy. 

The measured quantity of our system is the population 
imbalance between \g) and |e). For a given atomic state 
= a\g) + b |e), the population imbalance is defined as 
P = bb* — aa*. Keeping with the filter analogy, we can ex- 
press P as 



= ir 



dtr(t)A(t), 



(2) 



where P is the initial imbalance, A(i) = A + Su>(t), and 
r{t) is the impulse response [1], commonly referred to as the 
sensitivity function [16, 17]. The above equation defines r(t) 
and then calculating the effect of frequency noise on P reduces 
to calculating r(t). Inserting the above form of 5uj(t, to) into 
(2) and differentiating with respect to Aqb, we get 



OF 



dA(f> 



dtr(t)S(t-t ) = r{t ). (3) 
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FIG. 1. (Color online) Schematic diagrams, r(t), and |-R(/)| 2 for 
different spectroscopy sequences. An example total spectroscopy 
time of 60 ms is shown. Throughout this work, blue(green) repre- 
sents Rabi(Ramsey) spectroscopy and orange(red) represents l-(7-) 
echo pulse spectroscopy. The schematic diagrams shown in (a) sig- 
nify when f2(t) is nonzero by their "high" value and signify free 
evolution (£l(t) = 0) with their "low" value. For Rabi spectroscopy, 
the nonzero value of SI is selected such that the total pulse area of 
the spectroscopy pulse is n. For the other spectroscopy sequences, 
the nonzero value of Q is always 7r/(0.005) rad/s. The diagrams are 
offset in height for clarity. Plots of \R(f)\ in (c) and (e) are calcu- 
lated from the corresponding r(t) curves plotted in (b) and (d). Di- 
agrams and functions related to modified Ramsey and 7-echo pulse 
spectroscopy are plotted with dashed lines for clarity. 



We calculate P using (1) for a particular spectroscopy se- 
quence, r2(i), and with an added phase step of A</> occurring 
at time to- Finally, r(t) is calculated by allowing to to vary 
throughout the spectroscopy sequence. 

Eqn. (2) allows us to rigorously connect the observed fluc- 
tuations in P to A (7;). An essential quantity to consider is the 
variance of P, I 2 = (¥ 2 ^ — (P) , which can be expressed as 
[21]: 

p OO 

I 2 = (2tt) 2 df S„(f) \R(f)\ 2 . (4) 
Jo 

Here, R(f) is the Fourier transform of r(t) and S„(f) is 
the single sided frequency noise PSD of the laser in units 
of Hz 2 /Hz. Fig. 1 shows a schematic diagram for the spec- 
troscopy sequences used in this work along with their corre- 
sponding sensitivity functions, r(i) and |i?(/)| 2 , for a total 
spectroscopy time of 60 ms. 

For Rabi and Ramsey spectroscopy, the theoretical value for 
/ diverges since thermal noise, a fundamental limit to S v (/) at 
low frequency, has an / _1 character. For these measurements, 
we use the Allan Deviation [22] to characterize fluctuations in 
P. In particular, we consider the two-sample Allan variance, 



defined as 

7 ( 2 2) ^((P 4+1 -P 4 ) 2 ), (5 ) 

where the index, i, signifies the zth measurement of P. In the 
treatment of multiple measurements we consider the sensitiv- 
ity function as periodic with a period equal to the experimen- 
tal cycle time, T c . For this work, T c is approximately 1 s plus 
the spectroscopy time. I%s can also be expressed in terms of 
S v (f) in the following way [21]: 

/•OO 

I 2 = (2n) 2 / df ^(/)2sin 2 (7r/T c ) |i?(/)| 2 . (6) 
Jo 

Although the theoretical value of I( 2 ) remains finite for all 
data conditions, it does not properly account for coherent vi- 
brational or electronic noise that exists on our laser at frequen- 
cies above 20 Hz. This noise is aliased onto our measure- 
ments and leads to regular, slow oscillations of the measured 
P. Therefore, for echo pulse sequences we use / to charac- 
terize our measurements. The Rabi and Ramsey sequences do 
not suffer from aliasing of coherent noise since they do not 
have significant sensitivity to noise above 20 Hz for the spec- 
troscopy times used here. 

Our experimental setup follows that of our Sr clock [23, 
24]. Between 2000 and 3000 87 Sr atoms are cooled to about 
2 fiK in a one-dimensional optical lattice and nuclear spin po- 
larized into the ground 1 So rap — 9/2 state. The optical 
lattice is kept near the magic wavelength [25] for the \g) to 
|e) clock transition. Lattice-trapped atoms are excited with 
698 nm light according to the spectroscopy sequences shown 
in Fig. 1(a). The clock light propagates along the strongly 
confined axis of the lattice so that it probes the atoms in the 
well-resolved sideband regime, free from Doppler and recoil 
effects. Finally, the numbers of atoms in \g) and |e) are mea- 
sured to determine P. 

Using this setup, we have resolved 0.5 Hz spectral fea- 
tures [26, 27] and demonstrated the most stable optical clock 
[28]. These results are enabled by the ultra-stable laser that 
addresses the clock transition (hereafter termed "a-laser")- 
The stability of the a-laser is at its thermal noise limit of 
1 x 10 -16 fractional frequency units for ~1 to 1000 s. We can 
look for strong noise features at higher frequencies in the a- 
laser's S v (f) using an optical beat with a second laser (here- 
after termed "/3-laser"). The /3-laser has demonstrated thermal 
noise limited performance at the 10~ 15 fractional frequency 
level [4]. The optical beat (Fig. 2(a)) is limited by the thermal 
noise of the /3-laser out to 10 Hz, where it becomes limited 
by the noise floor of the detector; however, discrete features 
are visible above this noise floor. Pairs of narrow noise peaks 
are visible near both 22 and 30 Hz. Additionally, noise peaks 
are consistently measured at 24 and 60 Hz. The 60 Hz peak is 
dominated by detector noise and the 24 Hz peak is visible in 
previous beat measurements between two /3-lasers [4] . 

The spectroscopy sequences are performed such that the 
measured P is sensitive to frequency noise from the a-laser. 
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FIG. 2. (Color online) (a) The measured S u (f) of a beat between the a- and /3-lasers is plotted in red. A model S v (f) for the a-laser is 
plotted as a black band (see text), (b) Measured values of I or Im are plotted as a function of total spectroscopy time for the spectroscopy 
sequences shown in Fig. 1. Colored bands represent calculated values of / or 1(2) using the model S v (f) in (a) with QPN added in quadrature. 
A gray band indicates the mean calculated QPN for all experimental data. The black cross represents a measurement of the QPN (see text). 



For Rabi spectroscopy, we detune the a-laser from resonance 
by the half width at half maximum (HWHM) of the Rabi line 
shape and subject the atoms to a it -pulse. For Ramsey spec- 
troscopy, we tune the laser exactly on resonance and subject 
the atoms to two ir/2 pulses separated in time. We shift the 
phase of the final ir/2 pulse by ir/2 radians relative to the ini- 
tial 7r/2 pulse, which is equivalent to detuning by the HWHM 
of the central Ramsey fringe in the absence of phase shifts. 
The echo pulse sequences add to the Ramsey sequence a num- 
ber of 7r-pulses introduced such that the free evolution times 
between laser pulses are all equal. The echo pulse sequences 
act as a bandpass filter peaked at (n + l)/(2r) Hz, where n 
is the number of echo pulses and r is the total spectroscopy 
time. One can intuitively understand this behavior from the 
sensitivity functions in Fig. 1(d), which are periodic at this 
frequency. The Fourier transforms of r(t) explicitly demon- 
strate this frequency sensitivity in Fig. 1(e). To prevent the 
accumulation of pulse area errors, we switch the phase of the 
laser between adjacent echo pulses by ir rad. 

We use 80 consecutive measurements of the atomic pop- 
ulations to estimate the raw standard(pair) deviation and its 
statistical uncertainty. The standard(pair) deviation is divided 
by the measured contrast to give 1(1(2))- The contrast is de- 
termined by a fit to the measured line shape as a function of 
detuning for Rabi spectroscopy or a fit to the measured oscil- 
lations of population when the phase of the final spectroscopy 
pulse is scanned for Ramsey and echo pulse spectroscopy. 
Fig. 2(b) shows measured values of I or 7( 2 ) f° r different spec- 
troscopy sequences as a function of total spectroscopy time. 
Each data point represents a weighted mean of at least four 
measurements and the error bars are estimated from the vari- 
ance of the weighted mean. Different spectroscopy times are 
applied in a random order to avoid systematic drifts. Each data 
point consists of measurements separated by several hours to 
insure consistency of the data. 

We create a model S v (f) for the a-laser so that the calcu- 



lated I and 7(2) are consistent with experimentally measured 
values. The model S l/ (f), plotted in Fig. 2(a), accounts for a 
finite range of magnitudes for 1/f and white frequency noise 
which are necessary to simultaneously fit both the Rabi and 
Ramsey spectroscopy data. Due to its insensitivity to white 
and 1/f noise at the inferred magnitudes, the 7-echo pulse data 
is used to determine discrete noise contributions to the model 
S„(f). Two peaks are dominant in the 7-echo pulse data cen- 
tered at 0.135 and 0.180 s of total spectroscopy time corre- 
sponding to noise components near 30 and 22 Hz respectively, 
consistent with the optical beat experiment. The noise at 30 
Hz originates from alternating current motors and 22 Hz is 
approximately an acoustic resonance frequency of the lab. 

The functional form of the model S v (f) is 

h N ft- 
&(/) = fewhite + e ™ + ^ — ~ — ~2, (7) 



where hi, fi, and Tj are the height, frequency, and full width 
at half maximum (FWHM) for the ith noise resonance. The 
parameters for these resonances are given in [21]. We deter- 
mine /ithermai to be 1-5(4) x 10 -3 Hz 2 , consistent with the 
predicted thermal noise of the cavity [26], and ft w hite to be 
3.3(3) x 10~ 3 Hz 2 /Hz. The widths and frequencies of the 
Lorentzian functions included in the model S v (f) are chosen 
to be consistent with features measured in the optical beat and 
the magnitudes are determined by the measured deviation of 
P. A broad noise peak was also necessary in our model S u (f) 
to account for acoustic noise resonances between 40 and 80 
Hz. This noise creates a small broad increase in I between 60 
ms and 90 ms of total spectroscopy time for the 7-echo pulse 
data and was previously observed in a beat with the /3-laser 
before a large acoustic shielding box was installed around the 
optics and vacuum chamber of the a-laser. Although an exact 
relationship between S u (f) and a FWHM linewidth has been 
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FIG. 3. (Color online) (a) Measured S v (f) of an optical beat be- 
tween the a- and /3-lasers and (b) measured and calculated values of 
/ for a 7-echo pulse sequence with and without external modulation 
at 46.6 Hz. (c) Measured S v (/) of an optical beat between the a- 
and /3-lasers and (d) measured and calculated values of I for a 7- 
echo pulse sequence with and without feedback cancelation of noise 
resonances near 22 Hz. (d) also demonstrates the appearance of 60 
Hz noise due to a malfunctioning signal generator that increases the 
measured values of / at 0.06 and 0.2 s of total spectroscopy time. 



established [29], an analytic expression for the linewidth does 
not exist in the presence of 1/f frequency noise. Here, the 
observed linewidth depends on the measurement time [30]. 
By properly accounting for a finite measurement time, we nu- 
merically calculate the minimum observable linewidth of the 
a-laser to be 26(4) mHz at 698 nm, dominated by the 1/f com- 
ponent of our model S„(f) [21]. 

For each data point, the QPN is calculated for the mea- 
sured number of atoms and the mean excitation fraction. The 
mean QPN for each spectroscopy condition is calculated and 
added in quadrature with the calculated values of / and In) to 
more accurately represent the experimental data. These quan- 
tities are plotted as colored bands in Fig. 2(b) where the up- 
per(lower) edge of the band corresponds to the upper(lower) 
edge of the model S v (f). The mean and standard deviation 
of all calculated QPN values are represented in Fig. 2(b) as a 
gray band. The QPN is experimentally measured by the stan- 
dard deviation of repeated measurements of P following a 5 
ms, resonant, 7r/2 pulse. The measured and calculated QPN 
are consistent. 

To further test our working model we intentionally add 
noise to the a-laser. White noise is passed through a band- 
pass filter at 46.6 Hz with 2 Hz bandwidth and used to fre- 
quency modulate the a-laser with an acoustic optical modu- 
lator (AOM). Fig. 3(a) demonstrates the effect of the mod- 
ulation on the optical beat between the a- and /3-lasers. By 
adding noise at 46.6 and 93.2 Hz into our model S v (f), cor- 
responding to the 1 st and 2nd order contributions of the mod- 
ulation, we are able to fully account for the measured values 
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FIG. 4. (Color online) The solid blue (dashed green) line represents 
the 1 s stability limit for an optical clock due to the Dick effect using 
Rabi (Ramsey) spectroscopy as a function of total spectroscopy time. 
The blue dotted (green dot-dashed) line represents the QPN stability 
limit at 1 s for Rabi (Ramsey) spectroscopy assuming a collection of 
2000 uncorrelated atoms. Ramsey spectroscopy assumes 2.5 ms 7r/2 
pulses. The blue "x" denotes the implied single clock 1 s stability 
from [28]. 



of I with modulation. Calculated and measured values of / 
are shown in Fig. 3(b) for 7-echo pulse spectroscopy with and 
without 46.6 Hz modulation added to the a-laser. 

In addition to using the optical beat to confirm the validity 
of our atomic measurements, we are also able to harness the 
information contained within the beat to reduce discrete noise 
features on the a-laser. We filter the beat with a bandpass 
centered at 22 Hz with sub-hertz bandwidth. This signal is in- 
verted and fed back onto the a-laser with an AOM. We can ob- 
serve the effect of the 22 Hz modulation on the beat itself (Fig. 
3(c)). The effectiveness of this technique is demonstrated by 
a reduction in the measured / between 0.15 and 0.20 s of to- 
tal spectroscopy time when feedback is active (Fig. 3(d)). In 
order to reproduce the measured values of / with feedback, 
the magnitude of the 22 Hz features needed to be reduced by 
60% in the model S u (f) compared to the condition without 
feedback modulation. 

Having developed an accurate model for the a-laser' s 
S v (f), we can predict the stability this laser can achieve when 
used in an optical atomic clock. Here, the laser's frequency is 
slaved to the \g) to |e) transition in 87 Sr by periodic interroga- 
tion of this transition. The stability of this frequency standard 
is limited by the Dick effect [15-17], whereby the periodic 
nature of interrogation creates sensitivity to laser noise at har- 
monics of 1/T C . In Fig. 4, we plot the one second stability 
limit due to the Dick effect for Rabi and Ramsey spectroscopy 
as a function of total spectroscopy time, r, using the lower 
limit of the model S v (f). We assume a typical T c of 857.5 ms 
plus spectroscopy time. We find that for Rabi spectroscopy 
with 160 ms interrogation time, the Dick effect limits the sta- 
bility of a clock using the a-laser to 2.8 x 10~ 16 /y/r in frac- 
tional frequency units. For a comparison of two uncorrelated 
clocks, one operating with 1000 atoms and one operating with 
2000 atoms, we predict a stability of 4.6 x 10~ 16 jypr which 
is within 5% of the achieved stability reported in [28]. 

We thank R. Ozeri for stimulating discussions on dynamical 
decoupling and J. K. Thompson and Z. Chen for parallel work 
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DERIVATION OF THE DETUNING-DRESSED BASIS 
HAMILTONIAN 



We start with the Schrodinger equation for a two-level atom 
driven by a laser field with a frequency that can vary in time. 
This is most simply treated using the instantaneous phase of 
the laser, </>t(t) = ^u>L{t)dt, where Wi,(t) is the instan- 
taneous frequency of the laser. For a general state \if>) = 
c g (t) \g) + c e (t) |e), we get the coupled differential equations: 



get: 



ic g (t) 



>-c g {t) - ( e *iW + e -^W) Ce (t). 



(SI) 

Here, fi(t) is the Rabi frequency and us a is the frequency of 
the \g) to |e) transition. 



To simplify the solutions to these equations, we substitute 
the complex state coefficients with 

c e {t) ee e iu ^ 2 c e (t) and c g (t) = e - iuj " t/2 c g (t). (S2) 

Using this substitution, (SI) becomes 

<iL e -*w.*/ag B (t) + ie-**»*/%(t) 

= ^e^*/ 2 c e (t) - ^ (e«*) + e-'*W) e^*/ 2 c s (t), 

(S3) 

We can cancel an overall phase factor of e ~(+) iLJ a f / 2 from the 
top(bottom) equation in (S3) and simplify to get 



= _<M ( e <(#L(t)-«.<) +e -(%(t)+*'.t)j 



(S4) 



In order to proceed, we assume that the instantaneous laser 
frequency is the sum of a constant frequency and a small noise 
term that varies in time, u>L(t) = wl + 5uj(t). Then we can 
write the instantaneous phase of the laser as </>i(i) = u>L t + 
J* 8w(t')dt' ee uj Lo t + 8<j){t). Substituting this into (S4), we 



n(t) 

2 ' 

/ e i(wi, t+(5^(i)+a) a 4) _|_ e -i(wi o t+«0(t)-« a t)J g g (t), 



iCg(t) = \ 



= Q(t) . 
2 ' 

C (S5) 

Since, 50(4) is assumed to be small, we can make the usual ro- 
tating wave approximation where we neglect terms that rotate 
at the sum of ujl and oj a . The resulting equations are 



ice(t) = -^e-H At +**W^(t), 



(S6) 



where, A ee cjj, — w a- 

Another basis change is necessary before arriving at the fi- 
nal result, so we express (S6) in terms of new state coefficients 
defined as 



b e (t) 



-iAt/2~ 



c e (t) and b g (t) ee e tAt/2 5 g {t). (S7) 



After canceling overall phase factors we get the following 
differential equations. 



ib e (t) =-^e- i5 *(%(t) + fb e (t), 
ib g ( t ) =-m e iH(t) be{t) _A bg{t) ^ 

This gives that the Hamiltonian in this basis is 



(S8) 



H _ fl(t) ( e i5 ^ 
~H ~ ~~2~ ye-* 6 *® 

where a z is a pauli spin matrix. 



— 7r°~z 
2 



(S9) 



DERIVATION OF VARIANCE 

Here we derive the variance of P, starting with equation (3) 
in the main text. We assume that all measurements are made 
near P = so it will be sufficient to calculate the quantity 



(P 2 ) = 



OO /"OG 



dt 1 dt 2 A(* x )r(*i) A(*a)r(ta) 



oo J — oo 



(S10) 
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This expression can be simplified using the definition of the 
convolution operation along with the definition of the auto- 
correlation function, 

ft/(tl-t2) = </(tl)/*(*2)>, (SH) 

where Rf is the autocorrelation function of /. Then, equation 
(S10) simplifies to 

/oo 
dt{U A * r)(t)r(t), (S12) 
-oo 

where "★" represents convolution. Since an oscillator's noise 
properties are typically characterized in frequency space via a 
power spectral density, we employ Parseval's theorem, which 
relates the integral of two complex valued functions x(t) and 
y(t) of time to the integral of their Fourier transforms X(f) 
and Y(f) as follows 

/OO /"OO 
dtx(t)y*(t)= / dfX(f)Y*(f). (S13) 
-oo J —oo 

Here, /* denotes the complex conjugate of /. Since r(t) is real 
valued, we get that 

/OO pOO 
dt(R, A *r)(t)r(t)= / dfF[(K A *r)} (f)R'(f), 
-oo J — oo 

(S14) 

where T is the fourier transform operation and R(f) is the 
Fourier transform of r(t). The convolution theorem states 
that the Fourier transform of a convolution of two functions 
is equal to the product of the Fourier transforms of the two 
functions. Additionally, the Wiener- Khinchin theorem states 
that the power spectral density (PSD) of a function, /, is equal 
to the Fourier transform of the autocorrelation function of /. 
Thus, we get that 

/oo 
dfS A (f) \R(f)\ 2 , (S15) 
-oo 

where S A (f) is the double sided PSD of the laser in units 
of (rad/s) 2 /Hz. We can express the variance in terms of the 
more commonly used single sided PSD with units of Hz 2 /Hz 
to obtain the final result 

I 2 = (P 2 ) = (2^) 2 dfS v {f) \R(f)\ 2 . (S16) 
Jo 

Here, we have use the relationship, S A (f) = {2tt) 2 S v (/)/2 
for / > 0. 

DERIVATION OF THE TWO SAMPLE ALLAN VARIANCE 

We begin with the definition for the two-sample Allan vari- 
ance: 

I 2 2) = ^((P t+1 -P,) 2 ), (S17) 



Then the expression for the two-sample Allan variance is com- 
pletely analogous to (S10) under the substitution 

r(t)->r (2) (t) = r(t)-r(t-T c ). (S18) 

Since the Fourier transform operation is linear and a shift in 
time of T c simply adds a phase factor of e r 2lT ' lT ^S to the 
Fourier transform of r(t), we get that 

■F[r (2) (t)] =F[r{t)]-F[r{t-T c )] 

(S19) 

= 2ie-* iT °f sin (ttTc /)#(/). 
This shows that I 2 ^ can be expressed in terms of R(f) as 

POO 

1(2) = (2^) 2 / 2 sin 2 (tt/T c ) \R(f)\ 2 . (S20) 

Jo 



TABLE OF PARAMETERS FOR MODEL PSD FEATURES 



TABLE L S v (f) Resonant Features. 



Index 


fi (Hz) 


hi (Hz 2 /Hz) 


r 4 (Hz) 


1 


21.87 


1.2 [1] 


0.03 


2 


22.39 


0.6 [1] 


0.03 


3 


29.45 


0.15 


0.1 


4 


29.90 


0.08 


0.4 


5 


60 


0.012 
Added Noise 


27 


6 


46.6 


0.07 


2 


7 


93.2 


0.25 


2 



CALCULATION OF THE LASER LINEWIDTH 

To calculate the linewidth of the a-laser, we express the au- 
tocorrelation function for the laser's electric field, Re(t), in 
terms of the frequency noise PSD, S v (f), and take the Fourier 
transform of Re(t) to get the laser line shape function [2, 3]. 
We consider an electric field of the form 

E(t) = E e l{2 ™ L ° t+sm \ (S21) 

where vl = ljl /2it Then the Re(t) can be expressed in 
terms of £„(/) as follows: [3] 

R E ( T ) = E 2 e l2 ™ L ° T e~ 2j ~ S " (/) "° p' T) df . (S22) 

Since the model S v (f) that we deduce from our measure- 
ments contains a component proportional to 1//, the integral 
in the exponent diverges and we must take into account the 
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FIG. SI. (Color online) The numerically calculated FWHM of 
the a-laser as a function of observation time. Black circles repre- 
sent calculations using the model S v (f) that is quoted in the main 
text, excluding narrow resonance features that do not contribute to 
the FWHM. Blue diamonds(purple squares) represent calculations 
using the upper(lower) limit quoted for S u (f). The minimum value 
for each condition is outlined in red. These values determine the 
minimum observable linewidth for the a-laser: 26(4) mHz. 
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FIG. S2. (Color online) The calculated line shape, Se{Sv), using 
hwhito = 0.0033 Hz 2 /Hz, h th crmai = 0.0015 Hz 2 , and T = 44 
s. This measurement time corresponds to the minimum FWHM of 
Se{v) for these values of S u (f). A dashed gray line signifies the 
half maximum to guide the eye. 



flat rectangle function between t = and t = T . Therefore, 
the line shape observed over a finite time To will depart from 
a simple Fourier transform of Re(t) and take the form 



S E {v) = 2 1 I 1 - £M e- i2 ^ T R E (r)dr, (S23) 



m 

-To \ T 

which can be simplified to 



S B (Su) = 4E$ / 1 - j-i oob(2it8ut) 



T> 



(S24) 



x e 



-2f~ TQ SAf) ' in2 £ fT) df 



dr. 



where 8v = v — v^ Q . 

To calculate the minimum observable FWHM, we com- 
pute the FWHM of Se{5v) as a function of To for the model 
S v (/) quoted in the main text and its upper and lower bounds. 
The three conditions lead to three minimum linewidths, ob- 
servable at three different measurement times (Fig. SI) and 
determine the minimum observable linewidth of the a-laser 
to be 26(4) mHz. The calculated line shape, Se{5v), using 
h whito = 0.0033 Hz 2 /Hz and h therma] = 0.0015 Hz 2 is plot- 
ted in Fig. S2 for the measurement time that give the mini- 
mum FWHM, To = 44 s. In practice we only include the 
white, 1/f, and broad acoustic resonance #5 contributions to 
S v {f) for this calculation because the narrow peaks increase 
computation time and we have confirmed that these peaks do 
not significantly affect the calculated FWHM (< 0.1% ef- 
fect for tested times). For T < 30 s, the FWHM is domi- 
nated by the Fourier limit of the triangular windowing func- 
tion (FWHM = 0.8859/T) for S v {f) = 0) with a small con- 
tribution from white noise (/i w hitc = 0.0033(3) Hz 2 /Hz leads 
to a Lorentzian linewidth of 10(1) mHz). Between T = 40 s 
and To — 60 s, the line shape can cross 50% of its peak value 
multiple times: a result of side bumps in the line shape arising 
from the windowing function and from the increasing contri- 
bution of 1/f noise as To increases. We always calculate the 
FWHM from the lowest frequency half maximum crossing. 



measurement time over which the laser is observed. A mea- 
surement lasting T) s is insensitive to frequency noise below 
1/Tq Hz. This dictates that the lower bound of the integral in 
Eqn. (S22) should be set to 1/Trj, thereby making the integral 
finite. The finite observation time also scales Re{t) by a tri- 
angle bump function since the measurement is windowed by a 
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features was necessary for the data plotted in figure 3(d), corre- 
sponding to daily variation in the acoustic environment. A further 
60% reduction was necessary for the condition with feedback ac- 
tive (see main text). 
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